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We consider a two- level quantum system (qubit) which is 
continuously measured by a detector. The information pro- 
vided by the detector is taken into account to describe the 
evolution during a particular realization of the measurement 
process. We discuss the Bayesian formalism for such "selec- 
tive" evolution of an individual qubit and apply it to several 
solid-state setups. In particular, we show how to suppress 
qubit decoherence using continuous measurement and a feed- 
back loop. 



I. INTRODUCTION 

Studies of two-level quantum systems have acquired 
recently a new meaning related to the use of this sim- 
ple quantum object as an elementary cell (qubit) of 
a quantum computer!!] This paper addresses the mea- 
surement of a qubit state, so it necessarily touches the 
long-standing and still somewhat controversial problem 
of quantum measurement pa which is known under the 
name of quantum state "collapse" . 

Having in mind a solid-state realization of qubit (for 
different proposals see, e.g. RefsJjTj) let us emphasize 
that a realistic detector has a noisy output signal, so 
the measurement of a qubit state should necessarily have 
finite duration in order to provide an acceptable signal- 
to-noise ifatkL. In this situation the "orthodox" collapse 
postulate^itj cannot be applied directly, since the mea- 
surement is not instantaneous. The necessity of a more 
general formalism is obvious, for example, in the case 
when the qubit "self-evolution" changes the quantum 
state considerably during a measurement process. Even 
if there is no self-evolution, one can wonder what happens 
with the qubit state after a partially completed measure- 
ment (when the signal-to-noise ratio is still on the order 
of unity). So, we need a formalism to describe the grad- 
ual qubit evolution, caused by the measurement process. 
As will be discussed later, the Schrodinger equation alone 
is not sufficient for the complete description of this evo- 
lution, and should be complemented by a slightly gener- 
alized collapse principle. 

Continuous quantum measurement was a subject of ex- 
tensive theoretical analysis during last two decades, and 
there are two main approaches to this problem. One 
approach is based on the. .theory of interaction with a 
dissipative environment .E3£j Taking the trace over the 
numerous degrees of freedom of the detector, it is possi- 



ble to obtain a gradual evolution of the measured system 
density matrix from the pure initial state to the inco- 
herent statistical .mixture, thus describing the measure- 
ment process-Ejuj Since the procedure implies an aver- 
aging over the ensemble, the final equations of this for- 
malism are deterministic and can be derived from the 
Schrodinger equation alone, without any notion of state 
collapse. The success of the theory in describing many 
solid-state experiments has supported an opinion com- 
mon nowadays that the collapse principle is a needless 
part of quantum mechanics. Because of the dominance 
of this approach (at least in the solid-state community) 
we will call it "conventional" . 

The other general approach, to. continuous quantum 
measurement (see, e.g., Refs.EZrEEl) explicitly or implic- 
itly uses the idea of the state collapse. Since quantum 
measurement is a fundamentally indeterministic process 
so that the exact measurement result is typically unpre- 
dictable, the approach describes the random evolution 
of the quantum state of the measured system. The im- 
portant advantage in comparison with the conventional 
approach is the absence of averaging over the total en- 
semble; hence, it is possible to describe the evolution of 
an individual quantum system during a particular real- 
ization of the measurement process. The evolution of 
the measured system obviously depends on a particu- 
lar measurement outcome; in other words, it is selected 
by (conditioned on) the measurement result. So, this 
approach is usually called the approach of selective or 
conditional quantum evolution. There is a rather broad 
variety of fcrmalisms and their interpretations within the 
approach. EZl E3 Depending on the details of the studied 
measurement setup and applied formalism, different au- 
thors discuss quantum trajectories, quantum state diffu- 
sion, stochastic evolution of the wavefunction, quantum 
jumps, stochastic Schrodinger equation, complex Hamil- 
tonian, method of restricted path integral, Bayesian for- 
malism, etc. (for comparison between several different 
ideas see, e.g. Refo). The approach of selective quan- 
tum evolution is relatively well developed in quantum 
optics; in contrast, it was introducedjnto the context of 
solid-state mesoscopics only recently^ 

In the present paper, we. continue the development of 
the Bayesian formalismE3C3~E3 for selective quantum evo- 
lution of a qubit due to continuous measurement. Several 
issues of the formalism derivation and interpretation are 
explained in more detail than in previous papers. A new 
way of derivation is presented for a special case of low- 
transparency quantum point contact (tunnel junction) as 
a detector. We also discuss equations (briefly mentioned 
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FIG. 1. Tunnel junction as a detector of the electron 
position in the double-dot which affects the barrier height. 
The current I(i) (detector output) reflects the evolution of 
the density matrix pij(t) of the measured two- level system 
(qubit). 



in RehEj) for the evolution of a qubit measured by single- 
electron transistor, which go beyonjtLthe approximation 
used for a nonideal detector in Ref.El Special attention 
is paid to a regime outside the "weakly responding" limit. 
Finally, we discuss the operation of a quantum feedback 
loop which can suppress the qubit decohcrence caused by 
interaction with the environment. 



II. EXAMPLES OF MEASUREMENT SETUP 

The total Hamiltonian H. of a qubit continuously mea- 
sured by a detector, 



H — Hqb + Hdet + 'Hint, 



(1) 



consists of terms describing the qubit, the detector, and 
their interaction. The qubit hamiltonian, 



Hqb = - (c\ci 



c\c 2 ) 



H(c\c 2 + cjci), 



(2) 



is characterized by the energy asymmetry e between two 
levels and the mixing (tunneling) strength H (we assume 
real H without loss of generality). The Hamiltonian (||) 
is written in the basis defined by the coupling with the 
detector. We will refer to mutually orthogonal states |1) 
and 1 2) as "localized" states in order to distinguish them 
from the "diagonal" basis consisting of the ground and 
excited states, which differ in energy by U£l — (AH 2 + 



A. Double-dot measured by tunnel junction 



Our study will be applicable to several different types 
of qubits and detectors. As the main example we con- 
sider a double quantum dot occupied by a single electron, 
the position of which is measured by a low-transparency 
tunnel junction nearby (see Fig. |lj). Basically following 
the model of Ref.El let us assume that the tunnel barrier 
height depends on the location of the electron in either 
dot 1 or 2; then the current through the tunnel junction 



(which is the detector output) is sensitive to the elec- 
tron location. In this case the detector and interaction 
Hamiltonians can be written as 

Hdet = Eiajai + ^ E r ala r + ^ T(a\ai + aja r ), 



Lr 



AT 

Hint = ^ -y (qci - c\c 2 )(alai + aja r ), 



(3) 



where both T and AT are real and their dependence on 
the states in electrodes (l,r) is neglected. If the elec- 
tron occupies dot 1, then the average current through 
the detector is h = 2tt(T + AT/2) 2 p lPr e 2 V/h (V is 
the voltage across the tunnel junction and pi <r are the 
densities of states in the electrodes) while if the mea- 
sured electron is in the dot 2, the average current is 
h = 2tt(T - AT/2) 2 p lPr e 2 V/h. 

The difference between the currents, 



A/ = h - h 



(4) 



determines the detector response to the electron position 
(notice the different sign in the definition of AI used in 
Ref.Ej). Because of the finite noise of the detector current 
I(t), the two states of the system cannot be distinguished 
instantaneously and the signal-to-noise ratio gradually 
improves with the increase of the measurement duration. 
Let us define the typical measurement time r m necessary 
to distinguish between two states as theJime for which 
the signal-to-noise ratio is close to unityO 



2(A/) 2 



(5) 



where Si and S 2 are the low-frequency spectral densi- 
ties of the detector noise for states |1) and |2). (As 
will be seen later, r rn also determines the timescale for 
selective evolution of the qubit state due to measure- 
ment.) For a low-transparency tunnel junction Si, 2 = 
2e7i,2 coth(/3eI//2), where (3 is the inverse temperature. 
At sufficiently small temperatures /3 _1 <gC eV (we assume 
zero temperature unless specially mentioned) the detec- 
tor shot noise is given by the Schottky formula, 



S U 2 = 2el u 



(G) 



To avoid an explicit account of the detector quantum 
noise we will consider only processes at frequencies to <C 
eV/h (in particular, we assume t^ 1 <§; eV/K). 

The major part of the paper will be devoted to the 
detector in the "weakly responding" regime when two 
states of the detector differ only a little (one can also call 
this regime "linear" , while the term "weak coupling" is 
reserved for a different meaning), in particular, 



|A/|«/ , 7 = (h+h)/2, 
\Si-S 2 \4ZSo, S = (S 1 + S 2 )/2, 

so the typical measurement time is 



(7) 

(8) 
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FIG. 2. Single-electron transistor (detector) measuring 
the charge state of the single-Cooper-pair box (qubit). 

T. m = 2S /(AI) 2 . (9) 

For a weakly responding detector the timescale c/Iq of 
individual electron passages through the detector is much 
shorter than r m , so the current can be considered con- 
tinuous on the measurement timescale. 



B. Double-dot and quantum point contact 

Besides the low-transparency tunnel junction as a de- 
tector, we can also consider a quantum point contact with 
arbitrary transparency T which depends on the electron 
position in the double dot. This setup in the context of 
continuous quantum measurement has been extensively 
studied both experimentallyErcS and theoretically.ciJcS 
In spite of a somewhat different mathematical descrip- 
tion (we will not write the Hamiltonian explicitly) this 
case is very close to the case above, which we prefer be- 
cause of its simplicity. The obvious, new feature is the 
different formula for the shot noisejij 

Si, 2 = 2e/ 1 , 2 (l-T 1 , 2 ), (10) 

where — Tx^e^V / ~k%. Notice that for the quan- 
tum point contact as a detector we make the condi- 
tion (Q) for weakly responding regime a little stronger, 
|AJ| -C (1 — 7i : 2)^i,2, so that both transmitted and re- 
flected currents can be considered continuous on the mea- 
surement timescale. 



C. Cooper-pair box and single-electron transistor 

Another interesting measurement setup (Fig. ||) intro- 
duced in Rcf.E9 in the context of a solid-state quantum 
computer, is a single-Cooper-pair box measured by a 
single-electron transistor (a somewhat similar setup has, 
been recently used for the experimental demonstrationuil 
of quantum oscillations in the time domain). The qubit 
in this case is represented by two charge states of a small- 
capacitance Josephson junction. The Josephson coupling 
provides the matrix element H in Eq. (Q) which is as- 
sumed to be much smaller than the single-electron charg- 
ing energy, so that only two charge states (adjusted by 
the gate voltage to be close in energy) are important. The 



capacitively coupled single-electron transistor is sensitive 
to the charge state of the Cooper-pair box and serves as 
the detector; the current I{t) through the transistor is 
the measurement output. 

One can find the detailed discussion of the Hamiltonian 
for this measurement setup in Ref.EHl. The qubit state af- 
fects the energy of the middle island of the single-electron 
transistor (Fig. ||), so the interaction is of "density- 
density" type: 

Hint = {c\c 2 - c\a) (^^a^cim - const ^ , (11) 

where the factor in large brackets is the number of ex- 
tra electrons on the transistor island. In the ^tthodox" 
regime of sequential single-electron tunnelingBEj in the 
transistor, the energy change AE affects the rates of tun- 
neling through the two tunnel junctions and thus affects 
the average current /. In the simplest case when the 
electrons can tunnel only in a strict alternating sequence 
with rates Tl aricLT^, the average currents I\ and I 2 can 
be calculated astHl 

Ii = eTL,iTR,i/(TL,i + T Rti ), (12) 

where i = 1,2 corresponds to the charge state of the 
qubit. Well outside the Coulomb blockade range the dif- 
ference between the rates is Tl,2 — = —AE/e 2 RL 
and r flj2 — r fl „i = AE/e 2 R,R, where Rl(r) > h/e 2 are 
the resistances of tunnel junctions. 

The measurement time to distinguish between states 
|1) and 1 2) for this setup is given by Eq. (|^), in which 
the spectral density of the single-electron transistor cur- 
rent can be calculated using equations of RefsxaiEj (the 
Schottky formula used for this purpose in RefE3 is valid 
only in a limiting case) . In the special case corresponding 
to Eq. (|l|) the shot noise is given by the formulaEj 

Si = 2el t (Tl ti + I^J/Fls + T^) 2 . (13) 



D. Two SQUIDs 

One more solid-state realization of continuous quan- 
tum measurement of a qubit can be done using two flux 
states of a SQUID as a qubit and another inductively cou- 
pled SQUID as a detector .E3 The corresponding Hamilto- 
nian and calculations of the SQUID noise can be found, 
e.g. in Refo. A minor difference in the formalism is 
related to the fact that the typical output signal from 
a SQUID is voltage instead of current in the examples 
above. 



III. RESULTS OF THE CONVENTIONAL 
APPROACH 

The goal of the present paper is the analysis of a selec- 
tive evolution of the qubit state due to continuous mea- 
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surement, taking into account the detector output I{t). 
However, before that ie.t pas-review the results of the con- 
ventional approachE3'E2TE§E3 to this problem which does 
not take into account the detector output. 

We describe the quantum state of a qubit by the den- 
sity matrix pij in the basis of localized states |1) and |2), 
so that pa (pn + P22 = 1) is the probability to find the 
system in the state "i" if an instantaneous measurement 
in this basis is performed, while pu (P21 = P12) charac- 
terizes the coherence; in particular, |pi2| 2 = P11P22 cos-, 
responds to a pure state. In the conventional approacht3 
the evolution of p^ is calculated using the Schrodinger 
equation for the combined system including the detector 
and then tracing out the detectoc-dpgrjeea-of freedom that 
leads to the following equations HeacSEJ 

pu = -P22 = -2 — Impi2, (14) 

£ H 

Pu = i t P12 + i -r- (fin - P22) - T d P12, (15) 
n n 

where the effect of continuous measurement is described 
by the ensemble decoherence rate r^. (Such equations in 
similar problems when the environment causes dephasing 
are known for many years - see, e.g. Refs.E3'E30.) 

For a double-dot measured by a tunnel junction_/Fig. 
|l|) the decoherence rate has been obtained in RefS 



rr 2 f 



2c 



(16) 



Comparing this equation with Eqs. ^ and (^J) one can 
easily notice that Td has a direct relation to the typical 
measurement time r m : 



(2r m )" 



(17) 



This relation obviously remains valid in the weakly re- 
sponding regime when the decoherence rate can be ex- 
pressed as 



T d = (A/) 2 /45 



(18) 



In the case of aJuute-transparency quantum point con- 
tact as a detectoiocfl the ensemble decoherence rate has 
been mainly studied in the weakly, re s p onding regime. 
The most important for us resulted e3c3e3 is that for 
symmetric coupling Eq. (|l§h is still valid, just the shot 
noise is now given by Eq. ( |To| ) instead of Eq. (||) (as men- 
tioned, the temperature is zero). In the asymmetric case, 
if the phase of transmitted and reflected electrons in the 
detector is sensitive to states |1) and |2), then there is 
an extra term in the equation, fox, the decoherence rate, 
so the decoherence is fastcrOC3cZl'C3 than given by Eq. 

The inequality Td > (2r m ) 1 has been also obtained in 
Ref.El for a single-electron transistor measuring a single- 
Cooper-pair box. The interaction Hamiltonian (|ll| ) al- 
lows us to relate the dephasing rate, Td — (AE) 2 S m /4h 2 , 



to the low-frequency spectral density S m of the fluctuat- 
ing number m of extra electrons on the transistor cen- 
tral islfpjjd. These fluctuations have been calculated in. 
Refs.E3£3 within the framework of the orthodox theory.E3 
In particular, assuming the weakly responding regime 



and the two-charge- 
©-(H) we obtain 



: dynamics corresponding to Eqs. 



r rf = 



(A£) 2 r L r, 

h 2 (T L 

(notice a different expression in Ref. 

nor" 



(19) 



In this case 



2T d T„ 



sr 2 , 



+ il) 



(Tl + T R ) 2 (T 2 L h/e 2 R R - T 2 R h/e 2 R L f 



(20) 



so for T L ~ T R this product is necessarily large, T d » 
(2r m ) _1 , since Rr : l ~> h/e 2 . However, for very dif- 
ferent tunnel rates the dephasing can be comparable to 
(2r m ) _1 . Assuming Tr^> Tl one can simplify Eq. j2tj| ) 
to 



2T d r m = 8(T L /T R ) 2 (R L e 2 /h) 2 



(21) 



Formally this expression becomes less than one if 
e 2 r £ i?L < KT R /^/&] however, in this case the significant 
cotunneling makes the orthodox approach invalid-and the 
quantum noise contribution becomes important .L3 In the 
cotunneling regime (well below the Coulomb blockade 
threshold) r d should be obviously comparable to (2r m ) _1 
because in this case essentially the barrier height (the en- 
ergy of the virtual state) is sensitive to a measured state, 
so the detecting principle becomes similar to the case of 
Fig. |l|. The inequality Td > (2r m ) _1 should remain valid 
in the cotunneling regime as well; this fact will be obvious 
from the Bayesian formalism. 

The quantum backaction of a SQUID in ther4inear- 
response approximation was calculated in RefH. It 
was shown that the total energy sensitivity of a SQUID 
(e/ey — e^y) 1 / 2 (which takes the backaction into account) 
is limited by h/2. Here ey is the "output" energy sensi- 
tivity [the output signal of a SQUID is V(t)], ej describes 
the intensity of backaction noise, and e/y characterizes 
their correlation. From the inequality e/ey > h 2 /4 we 
easily get an inequality for spectral densities: s/sy > 
h 2 (dV/d§) 2 , where dV/d§ describes the SQUID response 
to the flux $. For the two- SQUID measurement setup 
considered in the present paper the qubit dephasing due 
to backaction noise is Td = (A$) 2 s//4?i 2 where A$ is 
the measured flux difference between two qubit states. 
Using the inequality above for the product s/sy we ob=. 
tain a lower bound for the ensemble decoherence rate:C3 
Td > (AV r ) 2 /4sy = (2T m )~ 1 similar to all other setups 
discussed above. This lower bound can be achieved only 
when the SQUID sensitivity is quantum-limited. 

Notice that the main equations (^4|)-(|l^) of the con- 
ventional formalism do not depend on the detector out- 
put I(t), and so they cannot be used for the predic- 
tion of the detector current behavior [for generality we 
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again choose the current as a detector output signal 
even though for a SQUID it should be changed to V(t)]. 
An important step towards this goal has been made 
in Ref£3 for a tunnel junction as a detector (a simi- 
lar analysis for tha-single-electron transistor has been 
performed in Refs.E3i£3). Let us divide the density ma- 
trix pij into terms corresponding to different numbers n 
of electrons passed through the measuring tunnel junc- 
tion, = J2 n Pij ( on ly diagonal terms in n are consid- 
ered). Then the evolution of these terms is given by the 
equations^ 

1 -2|lmp^ 2 , (22) 
1 +2|lmp^ 2 , (23) 



Pll 


- e Pll 


e 


P22 




e 




• £ n 


H 


P12 


= ^? 2 + 






y/hh 

e 


P12 



h + h 

2e 



P12 



(24) 

while Eqs. (0)-© can be derived from Eqs. (H)-® 
after summation over n. 

Even though these equations couple the evolution of 
the system density matrix with the number of electrons 
passed through the detector, they cannot predict the be- 
havior of the current I(t) and do not allow the calcula- 
tion of for a given realization of I(t). Actually, this is 
quite expected since the conventional formalism describes 
the ensemble averaged evolution while the analysis of a 
particular measurement realization requires a formalism 
suitable for an individual quantum system. (The use of 



the conventiona 
recent attempts 



ojmalism was the reason why several 
l£3'El to analyze the detector current 
were not very successful.) The analysis of a particular 
realization of the measurement process can be performed 
using the Bayesian formalism discussed in the next sec- 
tion. 



IV. BAYESIAN FORMALISM 

In the Bayesianjjopnalism (the name originates from 
the Bayes formulaol£3 for probabilities) which was de- 
rived only for the weakly responding (linear) regime, the 
evolution of the qubit density matrix during a particular 
measurement process is described by the equationsEj 

H 2AI . .. . . . 

Pu = -P22 = -2 — Impi2 + P11P22 [I{t) - 7oJ, (25) 



So 



. e .H. . 

Pl2 = 1 T- Pl2 + 1 -T- (Pll ~ P22) 

n n 

AI 



(Pll ~ P22)-£- [I(t) ~ 111] Pl2 ~ ldp!2 
Do 



(26) 



(in Stratonovich interpretation - see below), which re- 
place Eqs. ©-(Hi) of the conventional formalism. Here 



Id 



r rf 



(My 

45 



>0, 



(27) 



is the decoherence rate due to the "pure environment" 
only (ideal continuous measurement does not lead to this 
decoherence) . One can see that 7^ = in the example of 
a tunnel junction as a detector, which thus can be called 
an ideal detector, r\ — 1, where 



n 



1 



2± 



1 



2T d T„ 



(28) 



A similar ideal situation occurs for a quantum point 
contact when Yd — (AI) 2 /ASq, and also for the two- 
SQUID setup when the sensitivity of the measuring 
SQUID is quantum-limited and the output and backac- 
tion noises are uncorrelated. The important prediction 
of the Bayesian formalism is that in such an ideal sit- 
uation (which is experimentally accessible), an initially 
pure state of the qubit remains pure during the evolu- 
tion; moreover, an initially mixed state can be gradually 
purified in the course of continuous measurement H3 

Eqs. (^)-(Q) allow us to calculate the evolution of 
Pij for a given measurement output I(t). In order to 
analyze the behavior of I(t), these equations should be 
complemented by the formula 

I(t)-I = ^-( Pll -p 22 )+at), (29) 

where £(t) is a zero-correlated ("white") random process 
with the same spectral density as the detector noise,E£l 
= Sq . The stochasticity of the detector current does 
not allow us to predict exactly the evolution of pij in 
each particular realization of the measurement process; 
however, the formalism describes the mutual dependence 
of the stochastic evolutions of Pij(t) and I(t) and thus 
allows us to make experimental predictions not accessible 
by the conventional approach. 

When Eq. (||) is substituted into Eqs. (||)-(||), we 
get a system of nonlinear stochastic differential equa- 
tions. The analysis of such equations requires special 
care, since their solutionj-depends on the accepted def- 
inition of the derivativeEll (this happens because the 
noi se in creases with the decrease of the timescale, and 
so £ 2 dt = const = S^/2 does not decrease with dt). 



In Eqs. (|25|)-([26|) we have used the symmetric defini- 
tion, p(t) = lim T ^ [p(i + t/2) - p(t - r/2)]/r. This is 
the so-called Stratonovich interpretation of the nonlinear 
stochastic equations. The main advantage of this inter- 
pretation is that all standard calculus—formulas [for ex- 
ample, (fg)' — f'g + fg'} remain valid,E3 so the intuition 
based on usual (nonstochastic) differential equations typ- 
ically works well (this is the reason why we prefer the 
Stratonovich interpretation). Its other advantage is the 
correct limit in the case when the white noise term is ap- 
proximated, by a properly converging sequence of smooth 
functions O 

However, for some purposes (e.g., for averaging over 
stochastic variables and for numerical simulations) it is 
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more convenient to use another definition of the deriva- 
tive: p(t) = ]im T ^o[p(t + t) — p(t)]/r. This is called 
the Ito interpretation and it is the most commonly used 
interpretation in mathematical literature on stochastic 
differential equations. There is a simple rule of transla- 
tion between the two interpretations for an arbitrary 
system of equations 



i<(t) = G i (x,i)+F i (x,i)f(t) 



(30) 



in Stratonovich interpretation, the corresponding Ito 
equation which has the same solution is 



St 



+ F,(x,tK(t) 



dx k 



F k (x,t) 



(31) 



where Xi(t) are the components of the vector x(t), Gi 
and Fi are arbitrary functions, and the constant S% is the 
spectral density of the white noise process £(t). Applying 
this transformation to Eqs. (p5[), (|2^), and ( p9| ) we get 
the following equations in Ito interpretation: 



Pn 



-P22 
e 



-2 — lmp 12 
n 



2AI 

Pii/022— (32) 



Pl2 = lr Pl2 

n 

-(Pi 



H ( 



P22) 



P22)-5-Pl2 4W 



7d- 



S'o 



(Alf 
45 



P12 , (33) 



while the current I(t) is still given by Eq. (|29|). Similar 
equations (in a different notation) have been obtained 
in Ref.EH for a symmetric two-level system measured by 
an ideal detector (e = 0, jd = 0). Notice that the Ito 
interpretation has been used in the majority of theories 
describing-selective evolution due to quantum measure- 
ment (seeEZrEj and references therein). 

Using the Ito interpretation it is easier to see that aver- 
aging of the evolution equations over the random process 
£(t) (i.e. averaging over different detector outputs) leads 
to the conventional equations (|lj)-(]l^). However, for the 
analysis of an individual realization of the evolution, Ito 
equations are typically less transparent for physical inter- 
pretation. For example, the term —pi 2 (AI) 2 /ASo in Eq. 
([33|) does not actually cause decoherence in an individual 
realization but just compensates the noise term propor- 
tional to £ 2 dt due to the Ito definition of the derivative, 
and so pi2{t) does not decrease exponentially in time if 
H =/= 0. Similarly, the fact that the measurement tries 
to localize the density matrix in one of two states is not 
clear from Eqs. (|32|)-(|33"|) while it is obvious from Eqs. 

To avoid confusion due to the difference between 
Stratonovich and Ito interpretations, it is helpful to write 
the exact solution of Eqs. (|2^)-(|2^) [which is also the so- 
lution of Eqs. (|32|) (0)] in the special case H = 0: 



Pn(t + r) _ p n (t) exp[-(/(r) - Ii) 2 t/S ] 



P22(t + r) p 22 (t) exp[-(/(r)-/ 2 ) 2 r/5 ] 



(34) 



Pl2(t + T) 



pi2(t) e 



ier /fi 



[pn(t + T)p 22 (t + r)f 2 [ Pll (t) P 22(t)] 1/2 



where 



J(r) = - 



t + T 



/(*') dt' 



(35) 



(36) 



is the detector current averaged over the time interval 
(t, t + t). These equations have clear physical meaning: 
Eq. (|34|) is just the Bayes formula (see next section) while 
Eq. ( p5| ) describes gradual decoherence due to the "pure 
environment" characterized by 7^. 

A useful tool for analysis of the measurement process 
is Monte-Carlo simulation of an individual process re- 
alization. For this purpose we can use Eqs. (34)-j35|) 
complemented by the simulation of evolution due to fi- 
nite H. Let us choose a sufficiently small timestep At 
(much smaller than h/H) and apply the following algo- 
rithm. First, for each timestep (t, t + At) we pick the 

averag ed current / = (At)- 1 f f +lxt I(t') dt' as a random 
number using the probability distribution 



(2tt j D) 1 /2 

P22W 
(27TD) 1 /; 



exp[ 



2D 

(I-h 
2D 



(37) 



where D = So/ 2 At. Then I is substituted into Eqs. 
@-@ to calculate p tj {t + At) from Plj (t). The last 
step of the procedure is the additional evolution during 
At due to finite H (rotation in the pn-pi 2 plane). Then 
the whole procedure is repeated for the next timestep At 
and so on. 

An alternative algorithm can be based directly on the 
Ito equations (|32|)-(f33"|) which are more natural for nu- 
merical simulations than the Stratonovich equations be- 
cause of the "forward-looking" definition of the deriva- 
tive. For sufficiently small At (now much smaller than all 
timescales S /(AI) 2 , h/H, h/e, and 7^) we first calcu- 
late the averaged pure noise, £ = (At)- 1 J £(t')dt', 
as a random number using the Gaussian distribution 



P(0 = (2ttD)-^ 2 exp[-(0 2 /2£>], 



(38) 



where again D = So/2At. Then this number is substi- 
tuted into Eq. @: 

Pll (t + At) = p u (t) - 2At(H/H)lmp 12 (t) 

+ Pii(t) P22 (t){2AI/S Q )lAt (39) 

and similarly into Eq. ( |33"1) . Then the updating procedure 
is repeated for the next step At and so on. The detector 
current can be calculated using Eq. (p9|). 

Both Monte-Carlo algorithms are equivalent; however, 
the first algorithm is better because it allows longer 
timesteps. The equivalence for small At can be proven 
analytically using a second-order series expansion of Eqs. 
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FIG. 3. Solid lines: gradual purification of the qubit den- 
sity matrix p(t) in the course of continuous measurement, 
starting from the completely incoherent state. Dashed lines 
show the evolution starting from localized states, assuming 
the same detector current. 



(|34[)-(p5[) and has also been checked numerically. Notice 
that for At < S Q /(AI) 2 , the current distribution @ is 
indistinguishable from the distribution P(£ + AI(pn — 
p 22 )/2) given by Eq. Q. 

A typical result of the Monte-Carlo simulation is shown 
in Fig. [| The solid lines show a particular realization of 
the evolution of p(t) (diagonal and nondiagonal elements 
of the density matrix) for a symmetric qubit, e — 0, mea- 
sured by a detector with coupling C = Ti(AI) 2 / SqH = 0.1 
and ideality factor r\ = 0.7. The real part of pi2{t) is 
not shown since its evolution is decoupled from pn{t) 
and lmpi2(t). The completely incoherent initial state 
is chosen, pii(0) = 0.5, pi2(0) = 0. Nevertheless, the 
measurement leads to the gradual onset of quantum co- 
herent oscillations. This happens because the measure- 
ment randomly tries to localize the qubit, while the finite 
H provides oscillations when the state becomes at least 
partially localized. The qubit state is gradually purified, 
eventually reaching a pure state if the detector is ideal. 
For a nonideal detector (Fig. ||) the state remains par- 
tially incoherent, which decreases the amplitude of the 
oscillations. 

The qubit gradually "forgets" its initial state during 
the evolution and the density matrix p(t) becomes de- 
termined mostly by the detector record. To illustrate 
this fact, the dashed lines in Fig . |^ show the qubit evolu- 
tion calculated by Eqs. (f25[)-(26) starting from two local- 
ized states and assuming that the detector current (not 
shown) is exactly the same as in the measurement re- 



alization corresponding to the solid lines. As expected, 
after the time comparable to r m the dashed lines become 
close to the solid lines. 

The tendency to qubit state localization due to mea- 
surement can be described quantitatively using the de- 
terministic part of Eqs. ( p5| ) and (p9|). However, because 
of the equation nonlinearity the typical localization time 
Ti cannot have a unique definition. If we define it via an 
exponential-growth factor exp(£/r;) for pn (^-evolution 
when pu is close to 1/2, then 



n = 2S* /(A/) 2 



(40) 



which exactly coincides with the definition of the typical 
measurement time r m . [If for the definition we choose the 
exponential-decrease factor exp(— t/ri) when the state is 
almost localized, then r; would be twice smaller.] 



V. DERIVATION BASED ON BAYES FORMULA 



« derivation of the 
3, which was based 



In this section we briefly review t 
Bayesian formalism presented in Ref. 
on the correspondence between classical and quantum 
measurements. 

In the classical case (H = 0, pn = 0) the measure- 
ment process can be described as an evolution of prob- 
abilities pu and P22 which reflect our knowledge about 
the system state. Then for arbitrary At (which can be 
comparable to r m ) the average current / obviously has 
the probability distribution given by Eq. (|37]). After the 
measurement during At the information about the sys- 
tem state has increased and the probabilities p\\ and p 2 2 
should be updated using the measurement result I and 
the Bayes formula (p4]), which completely [-describes the 
classical measurement. [The Bayes formulaEirEj says that 
the updated probability P*(A) of a hypothesis A given 
that event T has happened in an experiment, is equal 

to P(A)P(J r \A)/Y, B l p ( B ) p ( :F \ B )] where p ( A ) is the 
probability before the experiment, P(F\A) is the condi- 
tional probability of event T for hypothesis A, and the 
sum is over the complete set of mutually exclusive hy- 
potheses.] 

The next step is an important assumption: in the 
quantum case with H = the evolution of p\\ and p 2 2 
is still given by Eq. (|34|) because there is no principal 
possibility to distinguish between classical and quantum 
cases, performing only this kind of measurement. Even 
though this assumption is quite obvious, it is not derived 
formally but should rather be regarded as a consequence 
of the correspondence principle. In other words, this is 
the natural generalization of the collapse postulate to the 
case of incomplete measurement. 

The comparison with classical measurement cannot de- 
scribe the evolution of pi 2 ; however, there is an upper 



limit: \p\2\ < [pnp 22 ] 1 ' 2 - Surprisingly, this inequality is 
sufficient for the exact calculation of pi 2 (£) in the impor- 
tant special case of an ideal detector and H = 0. Aver- 
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aging this inequality over all possible detector outputs I 
using distribution (j37j) we get the inequality 

\pia(t + T)\ < [pn(t)p22(t)} 1/2 exp[-(AI) 2 T/4S }. (41) 

On the other hand, for such averaged dynamics Ecj. flU] ) 
actually reaches the upper bound [see Eqs. (0)-(^5|) and 
@] in the cases discussed in Section III (tunnel junc- 
tion, symmetric quantum point contact, or quantum- 
limited SQUID as a detector). This is possible only if 
in each realization of the measurement process the ini- 
tially pure density matrix Pij(t) stays pure all the time, 
\pi2(t)\ 2 = pu(t)p2i{t). This fact has been thej-niain 
point in the Bayesian formalism derivation in Ref.Eij. 

As the next step of the derivation, a mixed initial state 
has been taken into account (for H = and an ideal de- 
tector) using conservation of the "degree of purity" [Eq. 
( |35| ) with jd — 0] which directly follows from a statistical 
consideration. The qubit state evolution due to finite H 
has been simply added to the evolution due to measure- 
ment. Finally, the interaction with the extra environ- 
ment (which does not provide any measurement result) 
has been taken into account by introducing the decoher- 
ence rate jd- 

First-order series expansion of the corresponding equa- 
tions for pij (t + At) leads to differential equations (p5|)- 
(^6|). The reason why we get equations in Stratonovich 
interpretation is that the first-order expansion is neces- 
sarily based on the standard calculus rules which are 
valid only in this interpretation. Using a second-order 
expansion we can obtain differential equations both in 
Stratonovich and Ito interpretations, depending on the 
definition of the derivative. 



VI. ALTERNATIVE DERIVATION OF THE 
FORMALISM 

Let us discuss now an alternative way of deriving 
the Bayesian formalism, which is based on Eqs. (0)- 
( pi| ) of the conventional approach (a somewhat similar 
derivation of the .Bayesian formalism has been recently 
presented in Ref.Ej). Since these equations have been 
derivedE£l only for the tunnel junction as a detector, we 
limit ourselves to this case. 

Eqs. (|22|)-(p4]) describe the coupled evolution of the 
qubit density matrix /Oy and the number n of electrons 
passed through the detector, considering the "qubit plus 
detector" as a closed system. We need to make a small 
but very important step in order to describe an individ- 
ual measurement process: we need to construct an open 
system which outputs classical information to the out- 
side. For this purpose let us introduce the next stage 
of the measurement setup which will be called "pointer" 
(see Fig. [|). By definition, the pointer deals only with 
classical signals while quantum description is allowed for 
the detector. 



Pii n (t) 



n(t k ) 



qubit 




detector 




pointer 




< — > 
T 


< — > 
T 


> 

T 



quantum 
interaction 



collapse 



classical 
information 



FIG. 4. The pointer is introduced into the model to ex- 
tract the classical signal from the detector. 

Let us consider the following model. The pointer does 
not interact with the detector most of the time, however, 
at time moments t = tk (k = 1, 2, . . .) the pointer mea- 
sures (in simple orthodox way) the total number n of 
electrons passed through the detector. By our assump- 
tion the measured n should be a classical number, so after 
each measurement by the pointer the number = n(tk) 
is well defined. However, during the "free" evolution 
of the "qubit plus detector" between the measurements 
by pointer the number n(t) gets smeared according to 
Schrodinger equation, i.e. satisfy Eqs. (p2|)-(f24j). By in- 
troducing sufficiently frequent readout (collapse) into the 
model we get the ability to describe the time dependence 
of the detector current. Of course, many other collapse 
scenarios are possible, however, if we show that within 
some limits the measurement process does not depend 
on the choice of times tk, this is a good argument justi- 
fying the generality of the model. 

The collapse at t = tk can be described in the orthodox 
way.EjlLS The probability P(n) to measure n electrons 
passed through a detector is 



(42) 



The measurement by pointer picks some random num- 
ber nk according to distribution (|42]), however, after 
the measurement this number is already well defined 
and the density-, matrix should be immediately updated 
(collapsed):!^ 



Pij(tk +0) = <*„,„„ Pij(t k + 0), 
Pij(tk + 0) 



Pu fc (*fc-0) + P&(**-0)' 



(43) 
(44) 



where S n ^ nk is the Kronecker symbol. After that the evo- 
lution is described by Eqs. (p2])-(p4[) until the next col- 
lapse occurs at t = tk+i- 

The detector current in our model has a natural av- 
eraging during time period between tk-i and tk and 
can be calculated as Ik = eAnk/ Atk, where Arik = 
n(tk) — n(tk-i) and Atk = tk — tk-x- Since the detector 
output is intended to reflect the evolution of the mea- 
sured system, tk should be sufficiently frequent, in par- 
ticular Atk -C h/H. For a while let us completely neglect 
the terms proportional to H in Eqs. (p2])-(p4]) and dis- 
cuss their effect later. Then these equations can be solved 
exactly. For the initial condition p™ (0) = <5 n .oPy(0) the 
solution is 
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= eM-ht/e) Pll (0), 



p? 2 (*) 



n\ 

Whht/e) n 



exp(-I 2 t/e) p 2 2(0), 



exp(- 



h 



2e 



■t 



iet 



(45) 
(46) 

)pu(0). (47) 



Similar equations describe the evolution after kth mea- 
surement by the pointer, just t is shifted by tk and n is 
shifted by n k - Using Eqs. (f43|)-([4"4|) we derive the itera- 
tive equations for the qubit density matrix: 



pn(tk) 



Pii(tfc-iUi A " fc exp(-7i Atfc/e) 
x [p 11 (t fc _ 1 )/ 1 A " fc exp(-/ 1 Ai fc / e ) 



P2a(*fc-i) exp(-7 2 At fe /e) 



P22(*fc) = 1 - /Oll(*fc) , 

Pi2(ife 



(48) 
(49) 



Pn{tk)p22(tk] 



PV2\tk-lj ,. 

l P ll{tk-l)p22{tk-l) 

x exp(ie Atk/K) , 



1/2 



while the probability P(rik) to get n 
(At fc /e) A " fc 



(50) 
rifc at t = tk is 



P(n k ) = 



7 An * exp(-/i At fe / e ) pu(t fc _x) 



(An*)! 



+ J 2 A " fc exp(~/ 2 Atfc/e) P22(t k -i)] ■ 



(51) 



It is instructive to check that the averaging of Pij(tk) 
over the result of measurement at t = tk gives simple 
equations 



Pn(ifc) 



/Oll(tfc-l), p22(tk) = P22(ife-l), 

Pi2(*fe) = /Oi2(ifc-i)exp(ieAt fc /?i) 



x exp[-(V/i - y/l a yAt k /2e 



(52) 



(53) 



which aro— consistent with the conventional equations 

dHH)H 

One can easily see that Eq. ([48| ) can be interpreted as 
the Bayes formula, while Eq. (pfl) i s the conservation of 
the "degree of purity" , similar to the approach reviewed 
above. The complete equivalence between Eqs. (^8|)-(|5l|) 
and Eqs. (|3^)-(|37|) is achieved if |AJ| <C Iq and also the 
probing time At* is much longer than the typical time 
Io/e between individual electron passages in the detector 
(so that the current is essentially continuous). In this 
case the Poissonian distributions (§^)-(fl6|) obviously be- 
come Gaussian, and so the probability distributions for 
the current I = eAn k /At k given by Eq. ([37]) and Eq. 
( |5l| ) coincide. Similarly, Eqs. (|4^)-(|50|) for py evolution 
coincide-with Eqs. (P4)~(P5h applied to an ideal detector, 

i d = oM 

If the probing period is within the range e/Io <C At* <C 
eI /(AI) 2 , the evolution of pij is smooth and so Eqs. 
(Es)-(E9h can be written in a differential form which co- 
incides with Eqs. (p5|)-(p6|) of the Bayesian formalism 
with H = and = 0. The effect of finite H can be 
now taken into account by the addition of obvious terms 



into Eqs. (|25|)-([26j). However, this can be done only if 
Atfc <C H/fr because in the opposite case the terms of 
more than the first power in H should be added into 
Eqs. (|4jf)— (p0|) indicating a nontrivial interplay between 
two effects. 

So, we have shown that in the weakly responding case, 
AI -C Iq, Eqs. (p2|)-([2^) of the conventional approach 
complemented by a sufficiently frequent readout (col- 
lapse), e/I < At k < min[eio/(A/) 2 , h/H] lead to the 
equations of the Bayesian approach.. The decoherence 
rate 7^ is zero because the modelEij describes a tunnel 
junction which is an ideal detector. 



VII. EFFECT OF COLLAPSE DUE TO POINTER 

The simple model considered in the previous section 
allows us to analyze the effect of the repeated measure- 
ments by pointer on the qubit dynamics in more de- 
tail and beyond the approximations of the Bayesian ap- 
proach. First, it is important to notice that in this model 
the event of collapse at t = tk does not disturb the qubit 
measurement by the detector. More specifically, the col- 
lapse with unknown result n k is equivalent to the absence 
of the collapse. To prove this fact, Eqs. (^)-(Q) can be 
averaged with the distribution ( [42|) that results in unity 
operator. 

The absence of disturbance by pointer is because in the 
model there are no density matrix elements which couple 
detector states with different number of passed electrons. 
Physically, this is a consequence of the assumption of 
low detector barrier transparency and infinite number of 
electrons in the detector electrodes, so that the "attempt 
frequency" is much larger than any collapse frequency 
(for a quantum point contact the necessary condition for 
this assumption is the large resistance, R ^> Ti/e 2 ). In 
other words, this model is intrinsically Markovian and 
the detector is classical in a sense that the passage of in- 
dividual electrons through detector is essentially classical 
(not quantum) random process JI3 

The absence of the disturbance by collapse with un- 
known result does not mean, however, that we can for- 
get about the collapse and make it only "at the end 
of the day". Any readout from the detector necessar- 
ily changes the qubit state (or in other words, informs 
us about the change) and thus affects the qubit evolu- 
tion. In the limit of sufficiently often readout, At k <§C 
mm(e/Ii,e/l2,h/H, h/e), the evolution equations (p^)- 
( p4| ) and (^2|)-(|4J) simplify because at most one electron 
can pass through the detector between readouts. During 
the periods of time when no electrons are passed through 
the detector, the evolution is essentially described by Eqs. 
(p2|)-(p4|) with n = 0, while the frequent collapses just 
restore the density matrix normalization, leading to the 
continuous qubit evolution: 



P11 



P22 



-2 — Impi2 
n 



AI 



■ P11P22, 



(54) 
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FIG. 5. The lines show a gradual purification of the qubit 
density matrix p(t) in the regime of quantum jumps (frequent 
detector readout with one-electron accuracy). The dots, tri- 
angles, squares, and crosses correspond to finite readout pe- 
riods Atk/{h/H) = 0.5, 1, 2, and 3 respectively. 
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P22) + ^(Pll 



P 22 )Pl 2 - (55) 



However, at moments when one electron passes through 
the detector, the qubit state changes abruptly; this 
change is given by Eqs. (|4^)-(|50|) with Arik = 1 and 
At k -> 0: 



hp n (t-0) 



Pll{t + 0) /i Pll (t-0) + Wt-0)' 
fta(* + 0) = l-p u (t + 0), 



Pia(i + 0) = pi2(t-0) 



pii(i + 0) p 22 (< + 0) 



pn(i-O) p 22 (t-0) 



1/2 



(56) 
(57) 

(58) 



and can be obviously interpreted as the Bayesian up- 
date. Equations (|5^)-(p8|) pspespond to the framework 
of "quantum jump" model.E3'L3 

It is easy to see that initially pure qubit state remains 
pure under quantum jump evolution (|5^)-(|58|) and the 
density matrix is gradually purified if started from a 
mixed state. The lines in Fig. || show a particular re- 
alization of such evolution for I\/e — H/fi, I1/I2 = 3, 
and completely incoherent initial state, pn(0) = 0.5, 
Pi 2 (0) = 0. Each discontinuity of curves corresponds 
to the passage of an electron through the detector (the 
jumps of pi 2 are typically smaller than the jumps of pn). 
The matrix element pn always jumps up because I\ > I 2 
and so the electron passage indicates that the state |1) is 
somewhat more likely than state |2). The jumps are more 
pronounced when pn is closer to 0.5 because the jump 
amplitude is Apn = A7piip 22 /(/ipii + I 2 p 22 ) [see Eq. 




At k /(»/H) 

FIG. 6. The average qubit coherence factor 9 as a function 
of the readout period Atk for the measurement process shown 
in Fig. §. 



(pitf)]. The model allows us to consider finite ratio /1/I2 
in contrast to Eqs. (|2^)-(|2^) of the Bayesian approach. 
In the limit of weakly responding detector, |A/| <C Io, 
the amplitude of quantum jumps (|5^)-(|5^) is negligible 
and Eqs. (p5|)-(p6|) are restored (inr-jthis sense they de- 
scribe a "quantum diffusion" modeEa). Notice, however, 
that equations of the Bayesian approach are applicable 
to a broader class of detectors. 

Since the model (|2^) -(|24| ) describes the ideal detector, 
the qubit state in Fig. ^| eventually becomes completely 
pure. However, if the readout period Atk is not suf- 
ficiently small, the information about the moments of 
electron passage through the detector is partially lost 
that decreases our knowledge about the qubit state. In 
the formalism this leads to a partial decoherence of the 
qubit density matrix. The symbols in Fig. |B| (dots, tri- 
angles, squares, and crosses) represent the readout with 
several different periods for exactly the same realization 
of a measurement process as for the lines which represent 
very frequent readout. When the readout is still suffi- 
ciently frequent (dots), we can monitor the qubit evolu- 
tion with a good accuracy (dots almost coincide with the 
lines). However, with the increase of the readout period, 
pn becomes close to 0.5 and pi 2 becomes close to zero, 
indicating a strongly mixed state. Figure ^ shows the 
corresponding decrease of the average coherence factor 
9 = 1 — 4(pnp 22 — |pi 2 | 2 ) with increase of the readout 
period Atk (equal time between readouts is assumed). 
The averaging is done over the readout moments for suf- 
ficiently long realization of the measurement process. We 
also tried few other expressions which describe the den- 
sity matrix coherence, all of them show a similar depen- 
dence on Atk- Notice the vanishing coherence in Fig. ^ 
when the ratio between Atk and the quantum oscillation 
period ttH/H is close to an integer numbe&n(the regime 
of quantum nondemolition measurementaJcil). 

In the special case H = all the information about 
the qubit state is contained in the result of the last mea- 
surement by pointer. This fact can be easily proven by 
applying Eqs. ([48|)-(|5"o|) twice and checking that result- 
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ing qubit density matrix does not depend on the result 
n\ of the first measurement while the dependence on the 
second result is the same as in the case of only one 
last measurement. Similarly, the probability distribution 
Pin-i) [see Eq. (^T|)| averaged over the result ri\ of the first 
measurement exactly coincides with P(n2) in absence of 
the first measurement. 

It is interesting to discuss the generalization of the 
model to the case of a low-transparency tunnel junction 
with finite temperature of electrodes. Then each of the 
currents 1\ and I 2 can be decomposed into two currents 
flowing in opposite directions, 



If II- = 0eV, 



(59) 



where i = 1, 2, j3 is the inverse temperature, and V is the 
voltage across junction. In this case Eqs. (|2^)-(|2~i|) are 
replaced by the following equations: 

■n _ H + 1 1 n , fl_ n-1 , h_ 
Pll — e Pll + Pll + Pll 

-2^Imp? 2 , (60) 

j+ _|_ j— j+ j- 

P22 = ~ P22 + ~ P22 + ~ P22 

+ 2^Im^ 2 , (61) 

n n -i-n n +i — (n n - n n ) - ^ + ^ + ^ + ^ o™ 
Pl2 — 1 ^ P12 + 1 ^ IPll P22^ 2 e ' 12 



nil 



ri-1 
Pl2 



n+l 



■Pl2 



(62) 



If the readout period At^ is much shorter than 
min(e// l ± , %/H), the detector still does not decrease the 
qubit coherence in spite of the finite temperature. How- 
ever, if individual electron passages are not resolved, 
the information about the number of electrons passed 
in each direction is lost that leads to the qubit deco- 
herence. In the framework of Bayesian formalism in 
the case of quasicontinuous current, e/lf <C Atk <C 
min[e/i/(AJi) 2 , h/H], we can easily calculate the output 
current noise Sq = 2eIo coth(/3eV/2) and the ensemble 
decoherence rate Td p,coth(/3eV r /2)(A/) 2 /8e/o (see also 
the derivation in Ref.E3) . Thus calculated detector ideal- 
ity factor, 



j] = [taurine V/2)] 5 



(63) 



becomes significantly less than unity at temperatures 
p- 1 > eV. 



VIII. DETECTOR WITH CORRELATED 
OUTPUT AND BACK ACTION NOISES 

Let us assume again a weakly responding (linear) de- 
tector and consider the case when the output detector 



noise is correlated with the "backaction" noise which pro- 
vides the fluctuations e(t) of the qubit energy level differ- 
ence and thus leads to the qubit dephasing. For example, 
this is the typ ical situation for a single-electron transistor 
as a detector .ESe3 In this case the knowledge of the noisy 
detector output I(t) gives some information about the 
probable backaction noise "trajectory" e(t) which can be 
used to improve our knowledge of the qubit state. The 
compensation for the most probable trajectory e(t) leads 
to improved Bayesian evolution equations SB 

H 2AI . , . t / \ 

P11 = -P22 = -2 — Impi2 + P11P22 -5— [I(t) - I \, (64) 
a 

. e . H . 

P12 = 1 r P12 + 1 (P11 - P22 ) 

- (pll - P22)^ [/(*) - Io] P12 

+ iK[I(t) - {puh + P22I2)] P12 - ldPi2, (65) 

where K = (de / 'd(p)Si v / Soh characterizes the correlation 
between the noise of current / through the single-electron 
transistor and the noise of its central electrode potential 
ip (Si v is the mutual low- frequency spectral density) .Lil 
The term in square brackets after K in Eq. ( |65| ) is just 
the "pure output noise" from Eq. (p9[). The decoherence 
rate jd in Eq. j6q ) is now decreased because of partial 
recovery of the dephasing: 



Id 



r„ 



(A/) 2 

45 n 



(66) 



The term containing K in Eq. fl6q) is proportional to the 
average (p(t) for given I(t). Performing ensemble aver- 
aging of this term [essentially, considering noise ipjt) as 
uncorrelated with I(t)], we can reduce Eqs. (|64|)-(|65|) to 
Eqs. (p5|) - ([2^), while additional ensemble averaging over 
I(t) leads to the conventional equations (|l4])-(|T5|). 

The obvious inequality 7^ > (in the opposite case 
the condition |pi2| 2 < P11P22 would be violated) imposes 
a lower bound for the ensemble decoherence rate IV 



r d > 



(A/) 2 K 2 S 
AS 4 



(67) 



which is stronger than the inequality 2TdT m > 1 (see 
Section III). 

Inequality ( |67| ) can be also interpreted in terms of 
the energy sensitivity of a single-electron transistor. 
Let us define the output energy sensitivity as ej = 
(dI/dq)~ 2 So/2C where C is the total island capaci- 
tance and dl/dq is the response to the externally in- 
duced charge q. Similarly, let us characterize the back- 
action noise intensity by e v = CS v /2 and the cor- 
relation between two noises by the magnitude e/ v = 
(dI/dq)~ 1 Si l p/2. Since in absence of other decoherence 
sources T d — S l p(CAE/2eh) 2 , where AE is the en- 
ergy coupling between qubit and single-electron transis- 
tor (see Section III), and using also the reciprocity prop- 
erty Ag = CAE/e = de/dip, we can rewrite Eq. (J67~ 
as 
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far " e? v ) 1/2 > h/2, (68) 

similar to the result of RefB (see also Refs.E&0). 
When the limit h/2 is achieved, the decoherence rate 



Id 



(Ml 
45 



{h/2f 



(69) 



in equations (J64|) — (J65|) for the selective evolution of an 
individual qubit vanishes, 7^ = 0. In this sense the de- 
tector is ideal, 77 = 1, where 



- = 1 _ Id = h 2 (dl/dq) 2 + (Siy) 



r 



SqSu: 



(70) 



even though it can be a nonideal detector (77 < 1) by 
the previous definition, 77 = h 2 (dl/dq) 2 / SoS v . [Notice a 
simple relation, 



1 



V = V 



(h/2) 2 _ 



(71) 



in absence of correlation between noises of <p(t) and I(t), 
{Si v ) 2 <C SqS v .} 

A similar conclusion is also valid for other kinds of de- 
tectors: a quantum-limited total energy sensitivity h/2 is 
equivalent-to detector ideality, 77 — l—Besides the tunnel 
junctionj23 quantum point contact ,E3'c3 and SQUID,Eil 
the regime of ideal quantum detection is also achiev- 
able by superconducting single-electron transisitoiH and. 
normal single-electron transistor in cotunneling rnpde.Lj 
(The resonant-tunneling single-electron transisitoiO has 
ideality factor comparable, but not equal to unity.) 



IX. QUANTUM FEEDBACK LOOP 

The Bayesian formalism allows us to monitor the evo- 
lution of an individual qubit using weak continuous mea- 
surement, thus avoiding strong instantaneous perturba- 
tions. This information can be used to control the qubit 
parameters e and H in order to tune continuously the 
qubit state in such a way that the evolution follows the 
desiredptrajectory (a similar idea has been discussed in 
RefsB^'Efl). This is possible even in the presence of deco- 
herence due to the environment and so presents an op- 
portunity to suppress such decoherence. 

Continuous qubit purification using a quantum feed- 
back loopH can be useful for a quantum computer. All 
quantum algorithms require the supply of "fresh" qubits 
with well-defined initial states. This supply is not a triv- 
ial problem since a qubit left alone for some time deterio- 
rates due to interaction with the environment. The usual 
idea is to use the ground state which should be eventually 
reached and does not deteriorate. However, to speed up 
the qubit initialization we need to increase the coupling 
with environment, which should be avoided. Another 
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FIG. 7. Schematic of continuous qubit purification using 
a quantum feedback loop. 



possible idea is to perform a projective measurement, af- 
ter which the state becomes well-defined. However, in 
the realistic case the coupling with the detector is finite, 
which makes projective measurement impossible. So, a 
different idea is helpful: to tune the qubit continuously in 
order to overcome the dephasing due to the environment 
and so keep the qubit "fresh" . 

The schematic of such state purification is shown in 
Fig. 0. The qubit is continuously measured by a weakly 
coupled detector, and the detector signal is plugged into 
Eqs. (|6^)"(|65|) [or into Eqs. (|25|)-(p6|) in a simpler case] to 
monitor the evolution of the qubit density matrix /ay (7;) . 
This evolution is compared with the desired evolution, 
and the difference is used to generate the feedback signal 
which controls the qubit parameters H and e in order to 
reduce the difference with the desired qubit state. 

We have simulated a feedback loop designed to main- 
tain the perfect quantum oscillations of a symmetric 
qubit (e = 0), so that the desired evolution is p\\ = 
[1 + cos(fir;)]/2, pi2 = isin(m)/2 where fl = 2H/h. Let 
us assume an ideal detector, 77 = 1, so that the qubit 
decoherence rate 7<j in Eqs. (|25|)-(|26|) is due to the ex- 
tra environment. The ratio between the decoherence rate 
and the "measurement rate" (A/) 2 /4So is described by 
the factor d = AS j d /(AI) 2 . 

To imitate a realistic situation the current I(t) is av- 
eraged with a rectangular window of duration r a run- 
ning in time, before it is plugged into Eqs. (p5|)-(p6|). 
So, thus calculated density matrix p a (t) differs (a lit- 
tle) from the "true" density matrix p(t) which is si- 
multaneously simulated by the Monte-Carlo method de- 
scribed in Section IV. The feedback signal is propor- 
tional to the difference between the desired oscil- 
lation phase Q(t — r a /2) and the phase calculated as 
4>(t) = arctajiplmpfa^)/^?!^)-^ {*))]■ Here the time 
shift r a /2 partially compensates the detector signal delay 
due to averaging. The feedback signal is used to control 
the qubit tunnel barrier: Hft,{t) = H[l — Fx Acf>(t — 7^)] 
where F is the dimensionless strength of the feedback 
and Td is an additional time delay (r^ = is preferable 
but not achievable in a realistic situation). 

Figure ^ shows typical realizations of the qubit 's pn 
evolution for C = h{AI) 2 /S H = 1, e = 0, r„ = O.lTi/H, 
Td = Q.05fi/H, and several values of the feedback factor 
F = 0, 0.3, and 3. No extra environment is assumed, 
d = 0. The qubit evolution starts from a localized state: 
Pn(0) = 1, Pi2(0) = 0, and the desired evolution is shown 
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FIG. 8. Particular realizations of the qubit evolution for 
a quantum feedback loop with strength F — 3, 0.3, and 
(thin solid, dashed, and dotted lines, respectively). Thick 
line shows the desired evolution. No extra environment is 
present, d — 0. 




t/(ft/H) 

FIG. 9. Operation of the quantum feedback loop (partic- 
ular realization) for several decoherence rates due to extra 
environment, d — j# x 4S'o/(A/) 2 = 0.3, 1, and 3 (thin solid, 
dashed, and dotted lines, respectively). Thick solid line is 
the desired evolution. F = 3, h(AI) 2 /S H = 1, e = 0, 

r a = o.m/H, r d = omn/H. 



by the thick solid line. Without a feedback (F = 0) the 
phase of quantum oscillations randomly fluctuates (dif- 
fuses) in time. However, for sufficiently large F the feed- 
back "locks" the qubit evolution and makes it close to the 
desired one. Further increase of F decreases the differ- 
ence between the actual and desired evolution. When 
F is too strong, the feedback loop becomes unstable. 
Overall, the behavior of this quantum feedback loop is 
similar to the behavior of a traditional classical feedback 
loop. In particular, we have checked that the increase 
of the averaging time r a and/or delay time eventually 
leads to synchronization breakdown. A decrease of the 
detector coupling C decreases the evolution disturbance 
due to measurement and allows more accurate tuning of 
quantum oscillations; on the other hand, in this case the 
feedback control becomes weaker and slower. 

Qubit decoherence due to the presence of an extra en- 
vironment prevents complete purification of the quantum 
oscillations so that the average qubit coherence factor 9 
becomes less than 100%. However, if the qubit coupling 
with the detector is stronger than the coupling with its 
environment, d%l, the feedback loop still provides qubit 
evolution quite close to the desired one (see Fig. g). Most 
noticeably, the phase of quantum oscillations does not 
diffuse far from the desired value Qt. So, for example, the 
spectral density of these oscillations has a delta-function 
shape at frequency fl (with exponentially small width) 
in contrast to the maximum value of 4 for the peak-to- 
pedestal f&jtit] in the case of quantum oscillations without 
feedback .E3'e3 



X. DISCUSSION 

The Bayesian formalism discussed in this paper 
presents j—Lels— any formalism of selective quantum 
evolutionEj lEI) a controversy in interpretation. First 
of all, a natural question is how it is possible that the 
qubit density matrix evolution can be described simulta- 



neously by the conventional equations (|14|)~(]15|) and the 
Bayesian equations (|2^)-(|26 Wand even also by the im- 
proved Bayesian equations (34])-(|65|)]. Which equations 
are correct? The answer is: all are correct depending on 
the problem considered. 

If only the ensemble evolution is studied (for exam- 
ple, the ensemble of particles is measured, as in typical 
nuclear magnetic resonance experiments) then the con- 
ventional approach is completely sufficient. It is also pos- 
sible to use the Bayesian equations; however, they should 
be averaged over all possible measurement results, after 
which they coincide with the conventional equations. So, 
the selective approach does not have real advantages for 
the study of the averaged evolution^ Chesides a significant 
computational gain in some casesEZrt£l) . There is still no 
advantage even for the majority of experiments with in- 
dividual quantum systems (see examples in Section II) if 
the averaging is done over a number of repeated exper- 
iments, disregarding the results of individual measure- 
ments (more exactly, when not more than one number is 
recorded as a result of each run). 

The principal advantage of the selective evolution ap- 
proach arises for continuous measurement of an individ- 
ual quantum system when the continuous detector out- 
put I(t) is recorded (or at least two numbers are recorded 
in each run). In this case the selective approach gives 
the possibility to make experimental predictions, unac- 
cessible for the conventional approach. The proposals of 
such experiments with solid-jstate qubits have been dis- 
cussed, for .-example, in Ref.ES for a one-detector setup 
and in Ref.H for a two-detector setup (the latter exper- 
iment seems to be realizable at the present-day level of 
solid-state technology). 

In this case the density matrices calculated by the con- 
ventional and Bayesian equations are significantly differ- 
ent. However, they do not contradict each other but 
rather the Bayesian-calculated density matrix is more ac- 
curate than the conventional counterpart. For example, 
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there are no situations when two approaches predict dif- 
ferent pure states of the qubit - then it would be possi- 
ble to prove experimentally than one of the approaches 
is wrong. Instead, in a typical situation the conventional 
equations give a significantly mixed state (so, essentially 
no predictions are possible) while the Bayesian equations 
give a pure state (and so some predictions with 100 % 
certainty are possible). A similar relation holds between 
Bayesian equations (p5|)-(p6|) and the improved Bayesian 
equations (pi|)-j6q): the latter give a more accurate de- 
scription of qubit evolution and allow us to make more 
accurate predictions. 

The difference between density matrices calculated in 
different approaches can be easily understood if we treat 
density matrix not as a kind of "objective reality" but 
rather as our knowledge about the qubit state (in accor- 
dance with orthodox interpretation of quantum mechan- 
ics). Then it is obvious that since Bayesian equations 
take into account additional information [detector output 
I(t)], they provide us with a more accurate description 
of the qubit state than the conventional equations. 

Another controversial issue is the state collapse due to 
measurement (here it is more appropriate to mention the 
mathematicalpfprmulation by LudersEilcPj rather than by 
von Neumannta) . The conventional equations are derived 
without any notion of collapse while the derivation of 
Bayesian equations requires either implicit or explicit (as 
in the model of Section VI) use of the collapse postulate. 
Philosophically, the collapse postulate is almost trivial: 
when the result of the measurement becomes available, 
we know for sure that the state of the measured system 
has changed consistently with the measurement result 
(even though it is generally impossible to predict the re- 
sult with certainty). In spite of being trivial, this pos- 
tulate in the author's opinion cannot even in principle 
be derived dynamically by the deterministic Schrodingcr 
equation because of the intrinsic randomness of the mea- 
surement result. In other words, the measurement pro- 
cess cannot be described by the Schrodinger equation 
alone because this equation is designed for closed systems 
while a quantum object under measurement is always 
an open system (even including the detector), since the 
measurement information is output to the outside world. 
(The incompatibility between quantum mechanics and 
"macrorealism" has been discussed, e.g., in Ref.U.) 

Following the orthodox (Copenhagen) interpretation, 
we can regard collapse not as a real physical process but 
rather as a convenient formal tool to get correct exper- 
imental predictions. In the author's opinion this tool 
is still irreplac eabl e (if we leave aside the many-worlds 
interpretations! 76 !! 77 !) for the complete description of the 
quantum realm. (Of course, in many cases the collapse 
postulate is not necessary as, for example, for the descrip- 
tion of decoherence due to interaction with the environ- 
ment - this problem has been solved with great success 
by the conventional approach.) 

Bayesian equations predict several quite counterintu- 
itive results. For example, even for a qubit with an infi- 



nite barrier between localized states, H = 0, the continu- 
ous measurement by an ideal detector leads to a gradual 
"flow" of the wavefunction between the states (for an ini- 
tially coherent qubit). The interpretation of this effect is 
rather difficult if we treat the wavefunction as objective 
reality; in contrast, there is no problem with the orthodox 
interpretation. Most importantly, experimental observa- 
tion of such effects in solid-state qubits is coming into 
the reach of present-day technology. These experiments 
would be extremely important not only for better under- 
standing of the foundations of quantum theory, but could 
be also useful in the context of quantum computing. 

The author thanks S. A. Gurvitz, D. V. Averin, H.- 
S. Goan, M. H. Devoret, A. Maassen van den Brink, G. 
Schon, and Y. Nakamura for useful discussions. 
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